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An accurate prediction of spacecraft avionics single event effect (SEE) radiation 
susceptibility is key to ensuring a safe and reliable vehicle. This is particularly important for 
long-duration deep space missions for human exploration where there is little or no chance 
for a quick emergency return to Earth. Monte Carlo nuclear reaction and transport codes 
such as FLUKA 12 can be used to generate very accurate models of the expected in-flight 
radiation environment for SEE analyses. A major downside to using a Monte Carlo-based 
code is that the run times can be very long (on the order of days). A more popular choice for 
SEE 1 calculations is the CREME96 3 deterministic code, which offers significantly shorter run 
times (on the order of seconds). However, CREME96, though fast and easy to use, has not 
been updated in several years and underestimates secondary particle shower effects in 
spacecraft structural shielding mass. Another modeling option to consider is the deterministic 
code HZETRN 2010 4 , which includes updates to address secondary particle shower effects 
more accurately. This paper builds on previous work by Rojdev, et al. 5 to compare the use of 
HZETRN 2010 against CREME96 as a tool to verify spacecraft avionics system reliability in 
a space flight SEE environment. This paper will discuss modifications made to HZETRN 2010 
to improve its performance for calculating SEE rates and compare results with both in-flight 
SEE rates and other calculation methods. 
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kilometers 
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I. Introduction 

An accurate prediction of spacecraft avionics single event effect (SEE) radiation susceptibility is key to ensuring 
a safe and reliable vehicle. This is particularly important for long-duration deep space missions for human exploration 
where there is little or no chance for a quick emergency return to Earth. Monte Carlo nuclear reaction and transport 
codes such as FLUKA 12 can be used to generate very accurate models of the expected in-flight radiation environment 
for SEE analyses. A major downside to using a Monte Carlo based code, though, is that the run times can be very long 
(on the order of days). 

A more popular choice for SEE calculations is the CREME96 3 deterministic code which offers significantly shorter 
run times (on the order of seconds). Elowever, CREME96, though fast and easy to use, has not been updated in several 
years and underestimates secondary particle shower effects in spacecraft structural shielding mass. Another modeling 
option which we investigate here is the deterministic code EIZETRN 2010 4 . E1ZETRN was developed primarily to 
support human space flight radiological health requirements, but it should also be possible to also use HZETRN for 
avionics SEE analyses. One advantage EIZETRN 2010 has over CREME96 is that it includes a more accurate 
representation of secondary particle shower products. 

In a previous paper 5 , we showed that the LET spectra generated by EIZETRN 2010 are limited to a maximum LET 
value of 11,596 MeV-cm 2 /g, reducing the utility of EIZETRN 2010 as a tool for calculating SEU rates. To eliminate 
this limitation, a relatively minor change was made to the EIZETRN 2010 source code to allow the calculated LET 
spectra output by EIZETRN 2010 to be extended to higher LET values. 

In the following paper, LET spectra calculations produced with CREME96 and a modified EIZETRN 2010 version 
are compared at different aluminum shielding mass thicknesses for both the interplanetary GCR environment and the 
International Space Station low-Earth orbit environment. Finally, space flight SEU rates are compared with rate 
predictions made using both deterministic codes. 


II. Methods 

The CREME96 GCR model was used to generate both the EIZETRN and CREME96 LET spectra data for this 
assessment. This was done to ensure that transport calculation results for the same input environment would be 
compared between the two codes. Solar minimum conditions were selected, which lead to higher intensity GCR 
exposure due to the inverse relationship between solar activity and GCR activity. An altitude of 362.5 km and orbit 
inclination of 5 1 .6° was specified for the LEO environment comparison. AP8MIN trapped protons were also included 
for the LEO environment definition. For the EIZETRN calculations, the CREME96 generated spectra were 
interpolated against the EIZETRN 2010 energy grid. The LEO (362.5 km) and free space particle environments are 
plotted in Figures 1 and 2, respectively. Only a subset of the heavy ion elements present in the environment models 
are plotted to avoid cluttering the two figures. 
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Figure 1. CREME96 solar minimum differential flux environment at 362.5 km and 51.6° inclination. 



Energy (AMeV) 

Figure 2. CREME96 solar minimum GCR differential flux environment for GEO and free space. 

As discussed in a previous paper 5 , the current version of HZETRN 2010 only includes elements out to Nickel (Z 
= 28). Extending the element range to higher Z particles would involve integrating a new environment model into 
HZETRN 2010 and was deemed out of scope for this study. The standard version of HZETRN 2010 also only outputs 
a maximum LET of 1 1,692 MeV-cm 2 /g for the GCR environment as compared to 30,364 MeV-cm 2 /g for CREME96 
(for max Z = 28). Extending the output LET range maximum in HZETRN 2010 can be accomplished by increasing 
the stopping power variable (IIS) in the HZETRN 2010 subroutines LETDIF and LETINT for calculating LET. For 
this study, IIS was increased from the default value of 700 to 780, which increases the maximum output LET to 7352 
keV/pm, which is equivalent to 31,552 MeV-cm 2 /g in silicon. 31,552 MeV-cm 2 /g is very close to the maximum LET 
of 30,364 MeV-cm 2 /g output by CREME96 for a maximum Z = 28. 

For the particle transport computational runs in the modified HZETRN 2010 and CREME96, spacecraft shielding 
was simulated using uniform aluminum shield layers of the following thickness: 0.1, 0.5, 1.0, 5.0, 10.0, 20.0 and 50.0 
g/cm 2 . A silicon detector behind the aluminum shield layer was used to simulate the sensitive volume of an electronic 
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device. Integral LET spectra were output from each transport code. These LET spectra were then compared between 
the two transport codes. Further analysis was then completed with CREME96 using the LET spectra from both codes 
to generate SEE rates for comparison. Results for 0.1, 5.0, and 50.0 g/cm 2 are presented herein. 

III. Results 

Integral flux LET spectra for shielding values of 0.1, 5.0, and 50.0 g/cm 2 generated from the modified HZETRN 
2010 and CREME96 with maximum Z = 28 are compared in Figure 3 through 5 for the LEO orbit and Figure 6 through 
8 for the GEO GCR environment. 



Figure 3. CREME vs. modified HZETRN LEO Integral flux at 0.1 g/cm 2 thickness. 



Figure 4. CREME vs. modified HZETRN LEO Integral flux at 5.0 g/cm 2 thickness. 
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Figure 5. CREME vs. modified HZETRN LEO Integral flux at 50.0 g/cm 2 thickness. 
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Figure 7. CREME vs. modified HZETRN GEO Integral flux at 5.0 g/cm 2 thickness. 



Figure 8. CREME vs. modified HZETRN GEO Integral flux at 50.0 g/cm 2 thickness. 

Radiation characterization parameters for several devices flown either on ISS, GEO, and deep space missions are 
given in 
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Table l 6 . The table shows SEE sensitive volume dimensions and Weibull fit parameters for each device. 
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Table 1. Device Parameters. 


Device 

Env 

CREME96 RPP 
x,y,z (pm) 

Onset 

(MeV- cm 2 /g) 

Width 

(MeV-cnf/mg) 

Exponent 

Limiting XS 
(pm 2 ) 

IMS1601EPI 

ISS 

39.5,39.5,5.92 

2750 

140 

0.95 

1560 

ISS SMJ4 16400 4Mx4 DRAM 

ISS 

1.05, 1.05,2.1 

420 

0.8 

1.7 

1.1 

ISS KM44S32030T 128Mbit SDRAM 

ISS 

2.42, 2.42, 0.24 

13000 

30 

1 

5.859 

ISS KM44S32030T 128Mbit SDRAM 

ISS 

1.25, 1.25,0.125 

14000 

30 

1 

1.563 

ISS KM44S32030T 128Mbit SDRAM 

ISS 

0.43, 0.43, 0.043 

1950 

30 

1.9 

0.186 

V4 XQR4VFX60 - BRAM 

ISS 

1.87, 1.87,3.74 

200 

70 

0.724 

3.5 

V4 XQR4VFX60 - Config. Memory 

ISS 

5.1, 5.1, 10.2 

500 

400 

0.985 

26 

V5 LX330T - Config. Memory 

ISS 

3.36,3.36,6.72 

500 

30 

1.5 

11.3 

Thuraya DSP Mega gate ASIC 

GEO 

2.5, 2.5, 1.76 

2700 

20.6 

1.2 

6.3 

Mercury Messenger ASIC 

GEO 

2,2,2 

300 

60 

6 

4 

Cassini OKI Solid State Recorder 

GEO 

6.32,6.32,6.32 

500 

32 

3 

40 

SOHO SMJ44100 4Mxl 

GEO 

7.07, 7.07, 2 

700 

15 

2.7 

50 

SOHO MHS CP65656EV 32kx8 SRAM 

GEO 

7.75,7.75,2 

1900 

17 

1.2 

60 

ETS-V PD4464D-20 64k SRAM 

GEO 

19, 19, 10 

500 

15 

2.9 

375 


SEE rates were calculated for these devices for three different LET spectra variations shown in Table 2 using the 
CREME96 HUP heavy ion SEE rate calculator. 


Table 2. LET Spectra Variations 


Code 

Max Z 

Max LET 
(MeV- cm 2 /g) 

Min LET 
(MeV- cm 2 /g) 

Modified HZETRN 

28 

31,552 

1 

CREME96 

28 

30,364 

1 

CREME96 

92 

101,000 

1 


To calculate SEE rates with CREME96 HUP using the HZETRN 2010 LET spectra, the spectra were first modified 
to conform to the CREME96 LET energy grid and then imported into the CREME96 website. SEE rates for shielding 
values for 0.1, 5.0, and 50.0 g/cm 2 are given in Table 5 for LEO environment devices and Table 6 - Table 8 for GEO 
and free space environment devices. All shielding depths are aluminum equivalent. Note that for the 50 g/cm 2 
thickness (Table 5 and Table 8) there is no data for CREME96 (Z=92) because CREME96 was unable to calculate 
SEE rates for this shielding thickness value. Rates for devices without the number of bits listed may be assumed to be 
device rates. 
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Table 3. SEE Rate for LEO Devices at 0.1 g/cm 2 


Depth = 0.1 g/cm 2 

SEE/bit-day 

LEO Environment Devices 

HZETRN 
Max Z=28 

CREME96 
Max Z=28 

CREME96 
Max Z=92 

IMS1601EPI 

1.79E-06 

1.76E-06 

1.77E-06 

ISS TMS44400 1Mx4 DRAM 

9.45E-09 

1.14E-08 

1.15E-08 

ISS SMJ4 16400 4Mx4 DRAM 

1.48E-07 

1.31E-07 

1.31E-07 

ISS KM44S32030T 128Mbit SDRAM 

1.27E-09 

1.33E-09 

1.34E-09 

ISS KM44S32030T 128Mbit SDRAM 

2.77E-10 

2.97E-10 

2.98E-10 

ISS KM44S32030T 128Mbit SDRAM 

3.38E-10 

4.09E-10 

4.10E-10 

V4 XQR4VFX60 - BRAM 

2.86E-07 

4.61E-07 

4.67E-07 

V4 XQR4VFX60 - Config. Memory 

1.23E-08 

1.10E-08 

1.10E-08 

V5 LX330T - Config. Memory 

9.29E-09 

8.82E-09 

8.85E-09 


Table 4. SEE Rate for LEO Devices at 5.0 g/cm 2 


Depth = 5.0 g/cm 2 

SEE/bit-day 

LEO Environment Devices 

HZETRN 
Max Z=28 

CREME96 
Max Z=28 

CREME96 
Max Z=92 

IMS1601EPI 

1.06E-06 

9.24E-07 

9.27E-07 

ISS TMS44400 1Mx4 DRAM 

8.77E-09 

1.07E-08 

1.08E-08 

ISS SMJ4 16400 4Mx4 DRAM 

1.04E-07 

9.02E-08 

9.04E-08 

ISS KM44S32030T 128Mbit SDRAM 

8.15E-10 

8.73E-10 

8.79E-10 

ISS KM44S32030T 128Mbit SDRAM 

1.76E-10 

1.94E-10 

1.95E-10 

ISS KM44S32030T 128Mbit SDRAM 

1.75E-10 

1.55E-10 

1.55E-10 

V4 XQR4VFX60 - BRAM 

8.37E-08 

6.30E-08 

6.33E-08 

V4 XQR4VFX60 - Config. Memory 

8.65E-09 

7.67E-09 

7.70E-09 

V5 LX330T - Config. Memory 

6.33E-09 

5.99E-09 

6.01E-09 


Table 5. SEE Rate for LEO Devices at 50.0 g/cm 2 


Depth = 50.0 g/cm 2 

SEE/bit-day 

LEO Environment Devices 

HZETRN 
Max Z=28 

CREME96 
Max Z=28 

IMS1601EPI 

2.15E-07 

1.39E-07 

ISS TMS44400 1Mx4 DRAM 

9.83E-10 

1.51E-09 

ISS SMJ4 16400 4Mx4 DRAM 

2.05E-08 

1.30E-08 

ISS KM44S32030T 128Mbit SDRAM 

1.13E-10 

1.22E-10 

ISS KM44S32030T 128Mbit SDRAM 

2.30E-1 1 

2.73E-1 1 

ISS KM44S32030T 128Mbit SDRAM 

3.47E-11 

2.31E-1 1 

V4 XQR4VFX60 - BRAM 

2.15E-08 

1.05E-08 

V4 XQR4VFX60 - Config. Memory 

1.65E-09 

1.10E-09 

V5 LX330T - Config. Memory 

1.09E-09 

8.49E-10 
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Table 6. SEE Rate for GEO Devices at 0.1 g/cm 2 


Depth = 0.1 g/cm 2 

SEE/bit-day 

GEO Environment Devices 

HZETRN 
Max Z=28 

CREME96 
Max Z=28 

CREME96 
Max Z=92 

Thuraya DSP Mega gate ASIC 

1.02E-07 

1.13E-07 

1.13E-07 

Mercury Messenger ASIC 

2.37E-1 1 

3.56E-1 1 

3.65E-1 1 

Cassini OKI Solid State Recorder 

4.94E-08 

5.96E-08 

5.99E-08 

SOHO SMJ44100 4Mxl 

1.38E-06 

1.49E-06 

1.49E-06 

SOHO MHS CP65656EV 32kx8 SRAM 

3.44E-06 

3.58E-06 

3.58E-06 

ETS-V PD4464D-20 64k SRAM 

6.49E-06 

7.15E-06 

7.16E-06 


Table 7. SEE Rate for GEO Devices at 5.0 g/cm 2 


Depth = 5.0 g/cm 2 

SEE/bit-day 

GEO Environment Devices 

HZETRN 
Max Z=28 

CREME96 
Max Z=28 

CREME96 
Max Z=92 

Thuraya DSP Mega gate ASIC 

3.50E-08 

3.93E-08 

3.94E-08 

Mercury Messenger ASIC 

7.91E-12 

1.40E-1 1 

1.43E-1 1 

Cassini OKI Solid State Recorder 

1.73E-08 

2.26E-08 

2.27E-08 

SOHO SMJ44100 4Mxl 

5.48E-07 

5.82E-07 

5.83E-07 

SOHO MHS CP65656EV 32kx8 SRAM 

1.56E-06 

1.58E-06 

1.58E-06 

ETS-V PD4464D-20 64k SRAM 

2.40E-06 

2.67E-06 

2.67E-06 


Table 8. SEE Rate for GEO Devices at 50.0 g/cm 2 


Depth = 50.0 g/cm 2 

SEE/bit-day 

GEO Environment Devices 

HZETRN 
Max Z=28 

CREME96 
Max Z=28 

Thuraya DSP Mega gate ASIC 

1.45E-09 

1.94E-09 

Mercury Messenger ASIC 

2.26E-13 

5.15E-13 

Cassini OKI Solid State Recorder 

7.33E-10 

1.04E-09 

SOHO SMJ44100 4Mxl 

4.25E-08 

3.40E-08 

SOHO MHS CP65656EV 32kx8 SRAM 

1.67E-07 

1.06E-07 

ETS-V PD4464D-20 64k SRAM 

1.43E-07 

1.42E-07 


Figure 9 compares in-flight SEU rates with rates calculated using F1ZETRN 2010, FLUKA and CREME96 using 
the accompanying least squares performance metric equations. Tabulated rates for in-flight, FLUKA, CREME96 and 
the modified F1ZETRN2010 code may be found in Table 9. The modified FIZETRN 2010 rate calculations are in 
general agreement with both the FLUKA and CREME96 SEE rate calculation results and compare favorably with in- 
flight measurements of SEU rates. 
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Figure 9. Comparing FLUKA, CREME, and modified FIZETRN 2010 SEU rate estimates with in-flight SEU rates. 


American Institute of Aeronautics and Astronautics 


11 


Table 9. Flight SEE Rates Compared to HZETRN and CREME96 Based SEE Rates. 


Device 

Env 

Median 
Shielding Mass 
g/cm 2 (3) 

In-Flight 
SEU/bit day 

CREME96 
Predicted 
SEU/bit day (1) 

HZETRN 
Predicted 
SEU/bit day 
(2) 

IMS1601EPI 

ISS 

34 (20) 

3.10E-07 

4.86E-07 

5.46E-07 

ISS SMJ4 16400 4Mx4 DRAM 

ISS 

10 

3.20E-09 

8.55E-09 

6.74E-09 

ISS SMJ416400 4Mx4 DRAM (4) 

ISS 

40 (50) 

3.70E-09 

1.30E-08 

2.05E-08 

ISS KM44S32030T 128Mbit SDRAM 

ISS 

40 (20) 

3.30E-10 

4.42E-10 

3.85E-10 

V4 XQR4VFX60 - BRAM 

ISS 

0.8(1) 

4.20E-08 

8.24E-08 

1.13E-07 

V4 XQR4WX60 - Config. Memory 

ISS 

0.8(1) 

3.80E-09 

9.04E-09 

1.05E-08 

V5 LX330T - Config. Memory 

ISS 

0.8(1) 

7.60E-09 

6.95E-09 

7.63E-09 

Thuraya DSP Mega gate ASIC 

GEO 

0.7(1) 

5.30E-08 

7.55E-08 

6.91E-08 

Mercury Messenger ASIC 

GEO 

1 

8.60E-10 

3.04E-11 

1.74E-1 1 

Cassini OKI Solid State Recorder 

GEO 

3.4(1) 

5.80E-08 

4.53E-08 

3.54E-08 

SOHO SMJ44100 4Mxl 

GEO 

1 

5.90E-07 

1.04E-06 

9.84E-07 

SOHO MHS CP65656EV 32kx8 SRAM 

GEO 

1 

1.70E-07 

2.59E-06 

2.56E-06 

ETS-V PD4464D-20 64k SRAM 

GEO 

5.8 (5) 

1.70E-06 

2.67E-06 

2.40E-06 

Table Notes: 

1) CREME96 Predicted (Z=92, Max LET -100 MeV-cm2/mg) 

2) HZETRN Predicted (Z=28, Max LET -32 MeV-cm2/mg) 

3) First value is the estimated on-orbit shielding for the device. Number in parenthesis is the shielding value used to calculate 
the CREME96 and HZETRN rates (if different from the estimated flight shielding). 

4) for this device shielding combination, the CREME96 prediction is for Z=28, Max LET -32 MeV-cm2/mg. CREME96 
would not calculate an SEE rate for the full Z=92 spectra. 

Median shielding mass and in-flight SEU rates taken from ref. 1 . 


IV. Discussion 

Direct comparison of the LET spectra produced by CREME96 and the modified HZETRN 2010 (Figure 3 -Figure 
8) shows that the modified HZETRN 2010 code produces LET spectra similar to CREME96. As shielding thickness 
increases, HZETRN 2010 predicts higher fluxes in the low LET range less than 1000 MeV-cm2/g. 

Comparison of CREME96 SEE rate calculation results shown in Table 6-Table 8 for GCR with maximum Z=28 
with results for GCR with maximum Z=92 shows that inclusion of the higher Z GCR elements had little impact on 
the overall SEE rates of the devices. This was expected, given the very low flux of GCR particles with Z > 28. Note 
also that all of these devices assessed have low or moderately low LET thresholds. Incorporating the high Z component 
of the environment is not so important for devices like these. The high Z elements become important for devices with 
higher LET thresholds and it is especially important for devices that may be susceptible to high-threshold destructive 
failures, such as single event latch-up and single event burnout. 

Figure 9 and Table 9 show that the HZETRN 2010 SEE rate predictions are generally conservative as compared 
to in-flight observations. HZETRN 2010 SEE rate predictions for the devices selected are also essentially the same as 
CREME96 SEE rate predictions. The disparity between the predicted and observed SEE rate for the Mercury 
Messenger ASIC is most likely due to the presence of a Tantalum (Z=73) overburden layer in the device which is not 
accounted for in the predicted SEE rate calculations. Nonetheless, the performance of HZETRN 2010 in predicting 
in-flight SEU rates is comparable to FLUKA and CREME96 for the devices studied here. 

What are the advantages and disadvantages of using HZETRN 2010 to model the SEE environment? The main 
advantage is that HZETRN2010 allows the use of any input environment model for calculating particle flux and LET 
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spectra at depth. Therefore, you can use any environment model to determine SEE rates. HZETRN 2010 also includes 
a more accurate representation of secondary particle shower products. 


There are two primary disadvantages to using HZETRN 2010 for SEE rate calculations. The maximum LET output 
by the standard version of HZETRN 2010 is 1 1,692 MeV-cm 2 /g. The maximum LET output by the modified version 
of HZETRN 2010 used in this paper is 30,364 MeV-cm 2 /g. To make HZETRN 2010 fully acceptable for spacecraft 
avionics SEE rate calculations, the range of the LET output tables needs to be extended to 100,000 MeV-cm 2 /g to 
handle high LET threshold devices. Additional modifications to HZETRN 2010 are required in order to properly 
transport the heavier ion species (Z>28) needed to produce these higher LET values. The other disadvantage is that 
HZETRN 2010 does not include a SEE rate calculation module, so the user would have to find an alternate method 
for calculating SEE rates. 


V. Conclusion 

The results in this paper show that a modified HZETRN 2010 is a viable alternative for determining LET spectra 
for SEU rate calculations for devices with low to moderate LET thresholds. To make HZETRN 2010 fully acceptable 
for spacecraft avionics SEE rate calculations, the range of the LET output tables needs to be extended to 100,000 
MeV-cm 2 /g to handle high LET threshold devices. The HZETRN 2010 user would also be responsible for performing 
their own SEE rate calculations since HZETRN 2010 does not include an SEE rate calculation module. 
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